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Abstract 



We use a variational wave function to calculate the energy of the interac- 
tion between external charges in the compact Abelian gauge theory in 2+1 
dimensions with mixed action. Our variational wave functions preserve the 
compact gauge invariance of the theory both in the vacuum and in the charged 
sectors. We find that a good estimate of the interaction energy is obtained 
only when we allow more variational parameters in the charged sector than 
in the vacuum sector. These extra parameters are the profile of an induced 
electric field. We find that the theory has a two-phase structure: When the 
D \ charge-2 coupling is large and negative there is no mass gap in the theory and 



no confinement, while otherwise a mass gap is generated dynamically and the 
theory confines charges. The pure Wilson theory is in the confining phase. 

I. INTRODUCTION 

The Abelian gauge theory in three dimensions has long been considered the archetype 
of a confining field theory. Polyakov [|l| demonstrated charge confinement with a classic 
analysis of the instanton plasma; the hope that the result would carry over to the instantons 
of Yang-Mills theory in four dimensions proved, unfortunately, to be unfounded Q . In the 
hands of Banks, Myerson, and Kogut ||, the theory became an exemplar of lattice gauge 
theory, displaying a single, confining phase from the trivial limit of strong coupling down 
to the expected essential singularity at weak coupling.^ It was also a stepping stone to 
elucidation of the four- dimensional Abelian theory, which confines at strong coupling but 



1 The corresponding demonstration for the four-dimensional non-Abelian theory still depends on 
numerical simulation 101. 
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contains only weakly interacting photons at weak coupling: The difference between the two 
phases is due to a condensate of magnetic monopoles at strong couping, closely related to 
the instantons of the three-dimensional theory. Subsequently, the Hamiltonian version of 
the theory became a laboratory for development of variational methods [||], providing in 
this formalism as well a specimen of nonperturbative confinement persisting down to weak 
coupling. 

In recent years there has been a revival of interest in variational approaches to non- 
Abelian gauge theories. Several calculations based on different variational ansatze have 
appeared in the literature 0,[7J. The approach of ||, a gauge invariant generalization of the 
Gaussian (mean field) variational approximation, was applied to the Abelian theory in ||, 
where it is clear that the wave function is a generalization of that used in || long ago. In the 
Yang Mills theory ||, there emerges a nonperturbative infrared mass scale, which makes the 
behavior of low momentum gluon modes in the vacuum very different from that calculated 
in perturbation theory. The comparison of the ultraviolet properties of this variational wave 
function with weak coupling perturbation theory was extended in p|,|T0 . 



The results presented in |6|,[],[I(]] are limited to the wave function of the vacuum state. 
The most interesting phenomenon in non-Abelian gauge theories, however, is confinement of 
fundamental charges. In a Lorentz invariant theory one could, in principle, calculate instead 
the vacuum average of a spacelike Wilson loop. In the Hamiltonian framework, however, 
Lorentz invariance is not explicit and most variational approaches will preserve it (at best) 
only approximately. The relation in the variational framework between the vacuum average 
of the spacelike Wilson loop and the energy of external charges is poorly understood. It is 
therefore all the more interesting to set up a variational calculation for the charged sector. 

In some recent papers |TT1 , [T2|] it was suggested that the charged sector can be described 
by the vacuum wave function, modified only to satisfy the new Gauss' law constraint. It was 
conjectured that further modification of the wave function in the presence of the external 
charges is unnecessary. We will refer to this state as the minimally modified vacuum. (Its 
exact definition will be given in Section 3.) It is far from evident that this simple procedure 
gives a good estimate of the energy of external charges. It does of course give a variational 
upper bound, but the quality of this bound may be poor. External charges introduce via 
Gauss' law an electric flux that is not present in the vacuum. It may be necessary to 
introduce additional variational parameters to allow this flux to spread out in the way that 
is most favorable energetically. 

The Abelian gauge theory in 2+1 dimensions is an ideal theater in which to examine 
this issue. The authors of || calculated the energy of an external dipole in this theory with 
the prescription advocated in [|rT| , |T2"| . The result was that the energy of the dipole is indeed 
proportional to its length and the string tension was determined. The vacuum expectation 
value of the spacelike Wilson loop was calculated in the same approach in |§. There it was 
found that the Wilson loop obeys an area law with a string tension consistent with the result 
of Polyakov [|. 

Superficially the results of H and M appear consistent. A closer look, however, reveals 
important differences. First, the vacuum state of || incorporates a dynamically generated 
mass whereas the correlations of magnetic and electric fields in the vacuum state of H have 
power law decays. Second, the string tension calculated in || (from the energy of external 
charges) is parametrically smaller than the result of M (from the spacelike Wilson loop). 
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Finally, a calculation of the spacelike Wilson loop with the wave function of || gives no area 
law. In the present paper we explain these discrepancies while providing a more general 
discussion of the calculation of the energy in the presence of external charges. 

We also take this opportunity to extend the variational analysis to the theory with mixed 
action. We define the theory with the lattice Hamiltonian 

H= l -a 2 Y,E 2 ni (a cos ga 2 B n + ^ cos 2ga 2 B t ] (1.1) 

Here a is the lattice spacing and the sums are respectively over the links and plaquettes of 
the two-dimensional spatial lattice. We constrain the constants a and (3 to satisfy 

a + P = \ (1.2) 

so that in the weak coupling limit, upon formal expansion to lowest order in g 2 , the Hamil- 
tonian reduces to the standard free Hamiltonian of 2+1 dimensional electrodynamics. We 
work in this weakly coupled regime, g 2 a <C 1, but go beyond perturbation theory with a 
variational wave function. 

The mixed Hamiltonian provides a new parameter for the theory. It has long been be- 
lieved that there is no non-confining phase in three dimensions, but generalized couplings 



along the lines of (|1 . 1| ) have not been adequately explored. Here we discover that this sim- 
plest of confining theories yet contains a surprise: a second-order confinement-deconfinement 
phase transition in the parameter a. We present in Section 2 a variational calculation that 
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determines the best variational vacuum of the theory. For a < ^ the energy is minimized 
for a state with a nonvanishing "mass gap," meaning that the correlation functions of the 
electromagnetic field decay exponentially at large distances. We will refer to this phase as 
massive. (This vacuum state is the same as that discussed in ||.) For a > \ no such 
mass gap is generated, and the correlation functions at large distances have a power law 
decay. This phase will be called massless in the following. (This state is the same as that 
considered in ||.) Note that the Wilson action a = 1 is in the massive phase.0 

We also calculate the expectation values of Wilson loops in the two phases. In the 
massive phase we find an area law for Wilson loops of arbitrary charge, with the string 
tension proportional to the square of the charge. In the massless phase the string tension 
vanishes and there is no area law. 

In Section 3 we calculate the potential between two external charges. First we consider 
the minimally modified vacuum according to the recipe of [||,|11],|T!J . We calculate the profile 



of the electric field and the energy of a dipole configuration. We find that the results of 
this calculation are inconsistent with what we would infer from the behavior of the spacelike 
Wilson loop. In the massive phase, where we found an area law, the electric field does not 
form a flux tube but rather takes the form of an ordinary, unconfined dipole field. (This is 



2 Morris |^] has studied the three-dimensional Abelian theory with generalized action and found 
phase transitions to phases with non-zero magnetic field. We have fixed a + /3 = 1 to prevent this 
in the weak-coupling limit we study. Our new, non-confining phase is unrelated to Morris'. 



3 



so even when the distance between the charges is much larger than the dynamical mass in 
the vacuum.) In the massless phase, amusingly, the situation is reversed. Even though there 
is no area law, the field far from the charges is screened and decays faster than a dipole field 
(but still as a power). 

The calculation of the energy gives the following results, again inconsistent with the 
results of Section 2. In the massless phase we reproduce the result of 0, a linear confin- 
ing potential. In the massive phase the energy of the minimally modified state diverges 
logarithmically in the infrared. 

These results convince us that the minimally modified state provides a poor upper bound 
for the energy of the external dipole.0 

In Section 4 therefore we introduce additional variational parameters into the trial state. 
These parameters are natural in the framework of the Gaussian approximation and corre- 
spond to a classical electric field present in the variational state. Varying the energy with 
respect to this classical field, we recover the results expected on the basis of the calculation 
in Section 2. In the massive phase the best variational state has the structure of an electric 
flux tube. Its energy is infrared finite and proportional to the distance between the charges. 
The string tension is the same as that calculated in Section 2. In the massless phase the flux 
spreads out into a dipole field, and correspondingly the energy has no linear dependence on 
the dipole size. 

Finally in Section 5 we discuss our results and their implications for a possible calculation 
of the string tension in non-Abelian gauge theories. 

II. THE VARIATIONAL VACUUM 
A. The wave function 

For a weakly coupled theory we expect the vacuum wave function to be not too different 
from that of a free theory. We therefore use a gauge invariant projection of a Gaussian wave 
function, as discussed in detail in ||. 

We present the formalism here in lattice notation, corresponding to the Hamiltonian 
(|1.1|). (We suppress the lattice spacing a henceforth.) Define a vortex field that satisfies 



(V x A 1 



—5 al>0 , V-A v = 0, (2.1) 



where n' is a site of the dual lattice. This is the vector potential corresponding to a magnetic 
field that is zero everywhere except at n' = 0, where it takes the value We demand that 
our variational wave function if) [A] be invariant under shifts A — > A + A v , which is an 
expression of the periodicity of H under B — > B + We also demand gauge invariance 
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Ben-Menahem [14] noted that the massless wave functions of H gave correlations inconsistent 



with the Euclidean analysis ||TJ,H| , and proposed a massive propagator to correct the discrepancy. 
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of the wave function. Hence we define a shifted field, shifted by a gauge transformation (j) n 
and by a vortex distribution m n /, 



A? m) = A n - (V0)n " E A V (n - n') 

n' 

or, for short, 

A (4>,m) = A-V<p- A v - m . 
Then the gauge invariant and periodic trial wave function is 

^[ A ] = E / [ d &J exp 



{"Vl 

Under a gauge transformation 



(4>,m) 
si 



*p[A + VA] = V[A] 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



since A can be absorbed in a shift in 0. The simple rotational structure of Gij = SijG 
that appears in the variational wave function ( |2.4j ) is consistent with perturbation theory, 
as discussed in ||. In this paper we take G{x) to be a real function. We now proceed to 
calculation of expectation values in the trial state (|2.4|) and the minimization of the vacuum 
energy. 



B. Normalization integral 

The normalization of is 

Z = = E / \W\W\ [dA] e -iAW.»)G-U(*."*) e -|^'.™ , )G-iA(*'.».') (2 6) 

{m,m'} 

We shift A by V0' + A v • m', and absorb the shift into <p and m, giving 

Z = E / [<*»A] e -^(*.™)G-^(*.») e -iAG-U (2 7) 

{m} J 

Now we combine the exponents according to 

A (4>,m) G -i A (4>,m) + AG" 1 A = 2A W2 ' m/2) G- 1 A W2 ' m/2) + -S{<f>, m)G- 1 S(<j), m) (2.8) 

2 

where 

S = V(j) + A v -m. (2.9) 

The first term in S has zero curl while the second is divergenceless. Furthermore, G^ 1 is 
translation invariant and proportional to the unit matrix. Thus SG~ 1 S has no cross terms 
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between m and <fi. We now shift A by V0/2 + A v ■ m/2, and all the fields decouple. We 
have then 



where 



Z — Zj^Z^Zy 



Za = det ttG 



1 \ 1/2 

det47r— G 



[2.10) 



(2.11) 



(2.12) 



z v = ex P 

{«V> 



V ft, 

Z v is the vortex partition function, with the vortex-vortex interaction D given by 
D(t' - s') =g 2 Y, AV ( r ~ r ') • G_1 ( r " s)A y (s - s ') . 

r,s 

The explicit solution of ( j2.1| ) is 

4» = (|^„.,„) 



so 

D = -^CT 1 ^ = -4tt 2 ^G- 1 . 
We can split off the r' = s' terms in (|2.13| ) and write 

1 



= £ ex P 

{"V> 



m r rD(r' — s')m*> 



HZ™*' 



where we have defined the vortex fugacity 



z = e 4 f 



(2.13) 



(2-14) 



(2.15) 



(2-16) 



(2.17) 



(2.18) 



In the interest of clarity, we adopt henceforth a continuum notation, indicating the 
ultraviolet cutoff a only where necessary. Moreover, since the interesting physics comes 
from the infrared, lattice effects can be approximated by a momentum-space cutoff A = a -1 , 
which can simplify formulas further. 

The variational function G appears both in the vortex-vortex potential and in the vortex 
fugacity. We expect the UV behavior of G to be the same as in the free theory, viz., 
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GT\k) ~ k 



(2.19) 



so 



D(0) 



^ d 2 k 4tt 2 



(2tt) 2 fc 2 



G-\k) 



and thus 



2ttA 



2! ~ e 9 a 



In the weak coupling region we have z < 1, justifying a restriction to m 
and ( gTTp . 



(2.20) 
(2.21) 

(2.22) 
0, ±1 in (P13|) 



C. Expectation values 

We calculate correlation functions of m as in || via a duality transformation. We add 
an z J ■ m term to the exponent in ( |2.13|) and use the formula 



e T a im Dm = const I [d~x\ e~ 



9 X-D -Xgix-Ti 



to obtain 



Noting thatQ 



we have 



Z V = J [d X ] e-^ D -^ + 2 cos( Xn + J»)] 



cos(x + J) = (cos x) : cos(x + J) : = z : cos(x + J) : 



Z„ = J [d X ] e-^ D '^ J] [1 + 2^ : cos( X + J) :] 

~ y -D^exp — g 2 xD~ l x + 2^ / gPx : cos + J(x) 

in continuum notation. Correspondingly ||, 

(m(ar)m(y)) = 2< 7 2 £r 1 (x - y) - V (/T^x) ir 1 *^) 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



The propagator of x is easily calculated. To first order in z, the only contribution comes 
from the tadpole diagrams, which have already been subtracted in ( p. 26 ). Therefore 



4 The normal ordering is performed relative to the free theory defined by the quadratic action in 
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d 2 xe lkx ( X (x)x(0)) 



2g 2 D- 1 (k) + 2z 
D{k) D 2 {k) 



2g 2 2g* 

The correlation function of the vortex density is then 



+ 0{z 2 



K(k) = / d 2 xe lkx (m(x)m(O)) = 2z + 0(z 2 ) 



(2.28) 



^2.29) 



which in this approximation does not depend on momentum. The k dependence will appear 
in the higher orders in z. 

Now we are ready to calculate the expectation value of the Hamiltonian ( |1.1|) . The 
calculation of the electric part of the energy is identical to that in ||. Using the definition 
( [2.4| ) we calculate 



V 



E 2 d 2 x 



E 



a 2 



if, 



(2, 
d 2 k 



G-\k)- 2 4z\^-k- 2 G-\k) 



9' 



(2.30) 
(2.31) 



The magnetic part is easily calculated since it has an exponential form and therefore with 
our trial wave function leads to a simple Gaussian integral. We find 



exp 



Ion/" d 2 k 
—n q I - — -7: 

4 y J (2tt) 2 



k 2 G{k) 



(2.32) 



The second factor is due to the vortices, and is different from unity only for odd values of 
n. Using (2.26) we find easily that (e t7Tm ) = e~ 4z . Expanding to leading order in g 2 and z, 
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l{ I- a 
— I a cos gB H — cos 2g B 



(2.33) 



(We have dropped an additive constant.) Finally, the expression for the variational vacuum 
expectation value of the energy is 



V 



(H) 



d 2 k 



G-\k) + k 2 G{k) - ^-z ( k~ 2 G- 2 {k) - A-a) 
q 2 V it 2 J 



(2.34) 



D. Determination of the ground state 



The expression (|2.34| ) has to be minimized with respect to G. From equation ( 2.18 ) and 
O) we find 
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Sz 



1 



5G(k) Ag 2 

The variation of ( 2.34 ) givesf] 

4vr 4 



— k- l G-\k)z 



(2.35) 



e-G- 2 (k) = ^ r zk- 2 G- 2 (k) [ 
9 J 



d 2 p 



p- z G-\p) 



7T Z 



[2.36) 



Eq. ( |2.36 ) has the solution 



where 



An 4 



G-\k) 



d 2 k 



m 



(27T) 



fc 4 

k 2 + m 2 



A;- 2 G- 2 (A;) - -L 

7T 2 



(2.37) 



^2.38) 



The main contribution to the integral in the gap equation ( [2.38] ) comes from momenta 
k 2 ^> m 2 . For these momenta k 2 G~~ 2 (k) = 1. We therefore see that (|2.38|) has a nontrivial 
solution when a < \. Using equations ( |2.18| ), ( |2.20| ), and ( |2.37| ) we obtain 



Air 4 



m 



71 

— exp - 

9 V 9 



d 2 p 



d 2 k 



(2n) 2 ^/p 2 + m 2 J J (2tt) 2 
which for g 2 <C 1 can be simplified to [cf. ( |2.22|) 1 



k 2 



k 2 + m? 7r 2 



a 



(2.39) 



m 2 = An 2 



71- 



-4a) 



7T 



exp 



(2.40) 



m is the mass gap of the theory, in the sense that it is the inverse of the spatial correlation 
length. Calculating, for example, the propagator of magnetic field, we find 



e igB me -igBr, 



,i 9 B\ 



e 2 



V 2 G(m-n) 



and at large distances (neglecting power-like prefactors) 

d 2 k 



V 2 G(x) 



(2tt) 



.(fc 2 + m 2 )V2 e 



2\l/2^,2k-x —mx 



(2.41) 



(2.42) 



This dynamically generated mass is Polyakov's result |Q, and is missing from ||. 

For a > ^r, ( |2.38| ) has no real solution for m. In this case the energy is minimized at 
the endpoint of the variational parameter range, m 2 = 0. We conclude that in the Gaussian 

2 

variational approximation, mixed action QED3 has a phase transition at a c = \. The mass 



We have dropped a term — ^p-zk G (k) from the right hand side of ( 2.36 ) is smaller by a 
factor of jJj- than the term retained, assuming G ~ k~ 4 at large k. 
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gap vanishes continuously at the critical point while (dm/da) a - 7^ 0. The phase transition 
therefore appears to be second order.0 

The energy density ( |2.34| ) contains ultraviolet contributions. The vacuum energy density 
in (almost) any field theory is ultraviolet divergent; thus it should diverge for any variational 
ansatz and indeed ( |2.34| ) contains such a divergence. The crucial question for us is whether 
the leading ultraviolet divergent terms depend on the variational parameter. This would 
indicate that the parameter is introduced in a way that directly affects the UV properties 
of the state. In other words such a dependence would mean that one is dealing with a set of 
states for which even the UV part of energy is not minimized yet. Such a variation would 
have little to say about IR physics. 

Examining (|2.34j) we see that the leading UV divergent terms are independent of the 
variational mass m. The UV behaviour of G is determined (in both phases) independently 
of m and is consistent with free field theory in the high momentum limit. (This indeed 
should be the case since the 3d U(l) gauge theory is superrenormalizable and its UV physics 
is governed by a fixed point at zero coupling.) Nontrivial information about the IR physics 
resides not in the leading terms in the energy but in the subleading ones. This makes us 
confident that the ultraviolet and infrared physics have indeed been separated in a clean 
way in our variational calculation. 

E. Spacelike Wilson loops 

Let us see whether charges are confined in the two phases. To this end we calculate the 
vacuum expectation value of the spacelike Wilson loop, 

W t [C] = ^exp (ilg j£ A • dxj \ = ^exp {ilg J B dS^j \ , (2.43) 

where I is an integer and the integral is over the area S bounded by the loop C . We have 

Wi[C\ = (j[e llmn ^ Z A l J DA exp (—AG' 1 A + ilg J^B dS^j (2.44) 

The second factor is a Gaussian integral, which gives the factor 

W A = exp ^ j s d 2 x J s d 2 y V 2 G(x - y)) . (2.45) 

In the limit of large S the leading behavior of the exponent is 

12 72 

- -g 2 S lim k 2 G{k) = --g 2 mS . (2.46) 



6 The value of the critical coupling a c = ^ should not be taken too literally, since G is modified 
from its continuum form close to the boundaries of the Brillouin zone. It is however clear that the 
phase transition is present in any cutoff scheme, at least in the present approximation. 
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In the massive phase this gives an area law with the string tensionT] 



a = jg 2 m . (2.47) 

In the massless phase this factor does not produce an area law. 

The first factor in equation ( |2.44 ) is different from unity only for odd I. It can be easily 



calculated to be 

W v = (j[e il7Tmn ^ = J D X exp -g 2 X D~ 1 X + 2zJ d 2 x :cos (x(x) - a(x)) : 



(2.48) 



where a(x) vanishes for x outside the loop and is equal to rr for x inside the loop. Recall 
that, from ( |2.16| ), 



D(k) = A<K 2 k- 2 G-\k) = An 2 {k 2 + m 2 )~ 1/2 . (2.49) 

At weak coupling we expand around a classical minimum of the exponent. In the massive 
phase, a < the inverse propagator D^ 1 is nonzero at zero momentum and thus it dictates 
the leading order solution xi x ) = 0. For this solution we find 

W v = e~ 4zS . (2.50) 



This is a sub-leading correction to the string tension ( |2.47| ), since g 2 m ~ e ~ K ^ 9 and z oc 



m 2 <C g 2 m. 

In the massless phase the situation is different. Here the kinetic term of the dual field x 
vanishes at zero momentum, £) -1 (0) = 0. Any constant value of x is therefore compatible 
with the kinetic term. The classical configuration that minimizes the vortex action in (|2.48|) 
for large area (S ^> z^ 1 ) is x = outside the area enclosed by the loop and x = 71 we U 
inside the area. The field interpolates between these two values in a layer of thickness z~ l 
in the neighborhood of the loop^. The action for this configuration clearly is proportional 
to the perimeter. Thus in the massless phase the string tension vanishes altogether. 

The spacelike Wilson loop tells us that the massive phase is confining with the string 
tension related in the expected way to the dynamically generated scale, a oc g 2 m. The 
massless phase does not exhibit an area law for spacelike Wilson loops and we therefore 
expect that the static charges in this phase are not confined. 



7 As was shown in |l5|] the dependence on I in this formula is incorrect. The correct result is a oc I 
rather than a oc I 2 . For multiply charged Wilson loops the nonlinearities of the compact theory 
are important and the Gaussian variational ansatz may be inadequate. This point is peripheral to 
the present paper and hence we will not pursue it. 

8 In the massive phase this classical configuration possesses an enormous action ~ g 2 mS, propor- 
tional to the area, since the kinetic term D" 1 does not vanish at zero momentum. The leading 
contribution to the Wilson loop therefore comes from configurations with % = well inside the 
loop as stated above. 
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III. THE MINIMALLY MODIFIED STATE 



As we mentioned in the Introduction, in a Lorentz invariant theory the string tension 
extracted from the spacelike Wilson loop should be the same string tension that determines 
the strength of the linear potential between two static charges. Lorentz symmetry, however, 
is not maintained in a variational calculation. Moreover, the theory we are studying is 
defined with an explicit ultraviolet regulator that itself breaks Lorentz invariance. It is thus 
interesting to perform a separate calculation of the potential between external charges. 



A. The wave function 

In the charged sector the gauge invariance condition (|2.5|) is replaced by 

^[A + VA] =V4A]exp(^][> n A n ^ , (3.1) 

where p n is a fixed, integer background charge distribution. (We have in mind a well- 
separated dipole.) A wave function that satisfies ( |3.1| ) is 



^[ A ] = J [ rf< ^] exp 

{™n'} 



r,s 



exp \ig^2pn(pn) • (3.2) 



The shifted field in ( |3.2| ) is defined along the lines of ( |2.3| ). 

A (4>,m) = A _ V( p _ A s . m _ ( 3>3 ) 

but now it is necessary to define the vortex field A s in a singular gauge. Like A v , the field 
A s satisfies 

(V x A 5 ) = — 5 n , ;0 , (3.4) 

V Jn> g 

but while A v is divergenceless we take for A s the solution where A is non-zero on the links 
dual to a string extending from (the plaquette of the vortex) to x = +oo, 

A nx = 

A s = / f for n x >0, n y = 
ny \ otherwise y } 

Our reason for using A s rather than A v in ( p.2|) is one of locality. A dipole can be 
created in the vacuum fl2.4Q by the string operator exp [ig ^4( z )], which places sources at 
x and y. This operator creates a string of electric flux taking integer values along the string, 
the most local way of preserving Gauss' Law when creating a dipole. A shift of A by A s 
commutes with this string operator, but a shift by A v does not. Shifting by A v will create a 
non-local, transverse electric field with fractional flux in addition to the string. In this light 
it appears unavoidable that the introduction of dynamical charges will immediately lead to 
a nonlocal and non-Lorentz- invariant theory. For more on this point see ||. 
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Another way of looking at it is to note that the difference between A v and A s can 
be absorbed by a shift in the integration variable r by ~ r')m r /, where 6{r — r') is 
the angle that the vector r — r' makes with the x axis. Using A v would therefore lead to 



an extra phase factor exp iY, r ,r' Pr9( r — r')m r > in the integral in ( |3.2|) . Under a shift of 
vortex density m this function is not invariant, but rather acquires a phase proportional 
to the charge density. The shift of vortex density can be viewed as a kind of large gauge 
transformation j§| and therefore this wave function belongs to a sector of the Hilbert space 
with a position-dependent angle." It is hard to imagine how this sector can define a local 
theory, especially in the presence of dynamical charges. 

Since A s differs from A v by a gauge transformation, it can be used interchangeably with 
it in the vacuum wave function (2^). The distinction is only meaningful in the charged 
sector, where it leads to a new vortex-charge interaction [see ( |3.9| ) below]. 

Taking G to be the same function as was used in Section |lj in the vacuum sector, this is 
the minimally modified state used by 
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B. Normalization and expectation values 

The normalization factor for the wave function ( |3.2| ) is 

Z = Z A Z^\p\Z v \p\ 



with 



Za = det nG , 

z Ap\ = ( det 471 ex P [~9 2 P ■ ^ G ■ P 

1 



L ~» ' ,2. 1 



z v[p\ = ex p 

{' m n' } 



\ r' / 



The new ingredient in (|3.9| ) is a vortex-charge interaction, 



■im r i 



r — r )m r i . 



(3.6) 

(3.7) 
(3.8) 

(3.9) 



(3.10) 



The vortex interaction potential D is again given by ( [2.16 ). 

First let us calculate the expectation value of the electric field in this state. Straightfor- 
ward algebra gives 



(Ei) 



d 

'dli 



di 

g^p G 

V 2 g 



(m) 



(3.11) 



The first term on the right hand side is the longitudinal field that gives the Coulomb inter- 
action between the charges in non-compact electrodynamics. The second term, entirely due 
to the compactness of the theory, is purely transverse as it must be in order not to spoil 
Gauss' law. 
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The calculation of the energy follows the same lines as the analogous calculation discussed 
in some detail in the previous section. For the kinetic term [cf. ( p. 30 )1, 



E 2 d 2 x 



7T~ 



+ — id xd yd G (x — y) (m(x)m(y)} 



(3.12) 



while the potential term is 



-l(«cos,5 + i^cos2,5)\ = i/^^)-^(Re(e-)-l) . (3.13) 



Putting these together we obtain for the energy of the minimally modified state 
(H) 



V r d 2 k r - ~2 

' \G-\k) + k 2 G{k) 



(2tt) 2 

.2 



g_ 

2 



1 

V 2 ' 



d 2 xp^-p 



7T 



2g 2 



dxdyd G (x — y) (m(x)m(y)} 



-- 2 Jd 2 x (Re(e i ™^)-l) . 



(3.14) 



To calculate correlation functions of the vorticity m, we again introduce the dual field x 
as in ( |2.23|) and obtain the relations 



(m) = -2iq 2 D- x ( X ) 
(m(x)m(y)) = 2g 2 D- 1 - Ag 4 (D~ l x(x)D- l x(y) 



where the averages of the dual field x are calculated with the partition function 
Z v = J Dxexp -g 2 xD~ l x + 2z J ^ x :cos (x(») ~ 6 



(3.15) 



(3.16) 



[The notation 9p represents the convolution of the source distribution p(r) with 6(r — r') 
In the weak coupling (small z) regime we solve the classical equation of motion, 



g 2 D- 1 X + zsm(x-9p) = 
We consider the massive (a < a c ) and massless (a > a c ) cases separately. 



(3.17) 



C. Field profile and energy in the massive phase 

First let us take a < a c . In this case the variational propagator G is massive. This mass 
is large compared to z since z oc m 2 <C g 2 m. Recall that D^ 1 ~ \/k 2 + m 2 . The classical 
equation ( |3.17| ) can then be solved in perturbation theory in To first order in z we find 

zD 

X = —sm(9p) (3.18) 
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and thus 



(m(x)} = — 2iz sin(9p(x)) 
(m(x)m(y)) connectcd = 2z cos(9p(x))5 2 (x - y) (3.19) 
/ e i7rm(x) ) =exp[-4zcos(0p(x))] 

Our first observation is that with this average value of the vortex density the electric 
field far away from the charges is not screened, but rather decays the same way as in the 
non-compact theory. To see this let us rewrite the expression for the electric field in terms 
of the dual field x, 

(Ek) = g^p + iLeydj (x) ■ (3.20) 

In order for the electric field to be screened away from the charges, the first and second 
terms must cancel. For this to happen the dual field x must behave asymptotically as 

x{x) — ► 9p(x) as x — ► oo . (3.21) 

To find the large-distance behavior of the solution ( |3.18[ ) we take for D{k 2 ) its value at zero 
momentum, -D(O) = (47i~ 2 m) -1 . We then find 

X(x -> oo) 1— sm(ep(x)) (3.22) 

g l m 

which is much smaller than required by ( |3.21| ), so there is no cancellation of the long range 
part of the electric field. 

Technically it is clear why we get the estimate ( |3.22| ) rather than ( |3.21| ) for the asymptotic 
behavior of the classical solution. If we were to minimize only the potential term in the dual 
action, ( |3.21| ) would be the result. In the massive phase, however, the kinetic term \D~ l x 
is the leading term for configurations dominated by low momentum components and it 
determines the infrared behavior of the solution. To illustrate this let us estimate the dual 
field action ( |3.16|) for a field configuration that behaves as ( |3.21|) . We take p to be a charge 



dipole with separation L. The main contribution to the action comes from the kinetic term 
in the infrared region, 

S x ~ g 2 mL 2 logV , (3.23) 



where V is the total volume (area) of the system. For the function ( |3.18 ) on the other hand 



the contributions of the kinetic and potential terms are of the same order and turn out to 
be 

S x ~ zL 2 logV. (3.24) 



This is much smaller than (|3.23j ) and therefore the configuration ([3.18 ) dominates the path 
integral over x- 

The electric field of a pair of static charges thus does not decay exponentially in the 
minimally modified state, even though the vacuum does contain a dynamically generated 
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scale. This leads us to suspect that the the minimally modified state is far from being the 
ground state in this sector. A calculation of the expectation value of the energy confirms 
this suspicion. The crucial terms are the last two terms, the vortex terms, in ( p.!4|) . To 
first order in z the contribution of this term to the energy of the minimally modified state 
relative to the vacuum is 

AE = m 2 g 2 J d 2 x [1 - cos(flp)] (3.25) 

Far from the charges, we have 8p(x) — > Lx\/x 2 when the dipole is parallel to the x\ axis. 
The integral in (|3.25|) is positive and logarithmically divergent in the volume. It is hard to 



believe that this infrared divergent energy has anything to do with the actual interaction 
energy of the charges. 



D. The massless phase 

In the massless phase, the kinetic term of the dual field vanishes at zero momentum, 
D~ l (k 2 ) ~ \k\. The solution therefore can not be expanded in powers of zD as before. The 
infrared behavior of the solution, however, can still be easily understood. It is determined 
by the potential term in ( 3.16|) and is given by ( p. 21 ). The distance at which this behavior 



sets in is determined by the value of momentum at which the kinetic and potential terms are 
comparable, x^ 1 oc \k\ = z/g 2 . The electric field is screened away from the dipole, decaying?] 
faster than 1/x. The "flux tube," carrying flux g in a region of transverse size g 2 /z, has 
energy 

AE ~ zL . (3.26) 

This is precisely the result of [H]. By itself it would suggest that the external charges are 
confined with the string tension a — z, despite the absence of an area law for the spacelike 
Wilson loop. We suspect, however, that just as in the massive phase, this contradiction 
indicates that the minimally modified state provides a very poor bound on the energy of the 
charges. 

A more reliable calculation of the interaction energy within the variational framework 
requires the introduction of additional variational parameters. This we do in the next section. 



IV. THE INTERACTION ENERGY IN THE EXTENDED VARIATIONAL 

ANSATZ 

As we have seen in the previous section the minimally modified state has a prescribed 
profile of the electric field. This profile, as it turns out, does not conform to one's intuitive 



9 The screening here is not exponential but rather power-wise. This is due to the nonanalyticity 
of the x propagator at zero momentum, z^O notwithstanding. Expanding D^ 1 around k = we 
find that the leading term D~ l ~ \k\, gives in coordinate space (xx) ~ 1/x 3 . 
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picture of the profile of the electric field in the charged state. When there is a dynamically 
generated mass and an area law for the spacelike Wilson loop, one expects the electric field 
between the external charges to form a flux tube. In the minimally modified state the field 
is just the Coulomb field of an electric dipole. In the massless phase, conversely, we expect 
to find a dipole field, but instead we obtain a screened distribution similar to a flux tube. 
Moreover, the energy of the minimally modified state is unnaturally high. 

A. The extended ansatz 

The above leads us to believe that the minimally modified state is very far in its physical 
properties from the minimal energy state in the charged sector. It is apparent that the main 
missing feature is freedom for the electric field to adjust its profile. We therefore introduce 
this additional variational freedom in our ansatz, writing 

^[A]= E /[^nlexpL^E^G-^r-^^HexpLEPn^ 

{m n ,} L A r,s J \ n / 

xexp^E e n" A«" m) ) • (4.1) 



We take the classical background field e to be transverse, V • e = 0, and we will treat it as 
an additional variational parameter in the charged sector. 

Retracing the calculations of the previous sections we find that the normalization Z is 
modified only in the vortex factor fl3.9|) , wherein the external field h r > acting on the vortices 
becomes 

27T 1 

h = ——e ij d i e j + e P . (4.2) 

The expectation value of the electric field is 

(E i ) = g^ s p-jG-^{m)+e i . (4.3) 

[Note that e enters ( |4.3| ) as a simple shift and also as part of the external field ( f4.2| ) that 
determines (m).} The energy of this state is 

+ 2(f J ^ X( ^y® 2( ~* 2 ( x ~v) ( m ( x ) m (v)) 

-j 2 Jd 2 x (Re(e l ™^)-l) 

+ l -QG- l (m)+ l -j<Pxe\ (4.4) 
where the penultimate term contains the potential C( x )i defined via 
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e * = h €ijdjC • (4 - 5) 

Thus we have also 

h = -C + 9p. (4.6) 
The duality transformation now results in the following Lagrangian for the dual field \i 

L = g\D-\ ~2zJ d 2 x : cos ( X (x) +(-6p): . (4.7) 

The correlation functions of m are related to the correlation functions of \ via ( |3.15| ). 

B. The massive phase 

In the massive phase, we obtain as before to first order in z 

(m(x)) — — 2izsin(9p(x) — ((x)) 
(m(x)m(y)} = 2z cos(9p(x) — ((x))5 2 (x — y) (4.8) 
'e mm{x) ) = exp[-4^cos(#p(x) - ({x))\ 



The e-dependent piece of the energy is 

AE = m 2 g 2 J d 2 x [1 - cos(9p(x) - C(x))] 



+z / d 2 xd 2 yC(x)G-\x - y) sin(9p(y) - ((y)) 

J d 2 x(V 2 ( (4.9) 



8tt 2 



The quantity ( (4.9| ) is to be minimized with respect to (. It is obvious immediately that we 
will have ( — > 9p at large distances from the sources, so that the infrared divergence found 
in the minimally modified state will disappear. Noting that the second term is formally of 
order g 2 relative to the first, we drop it for now, subject to a consistency check at the end. 
The minimization equation for ( then becomes very simple, 

V 2 C - m 2 sin(C - Op) = . (4.10) 



To study (|4.10 ) it is convenient to define ( = £ — 9p which satisfies the sine-Gordon 



equation with a singular source term, 

V 2 C-m 2 sinC = S (4.11) 

The source term S consists of a dipole layer along the line between the external point charges. 
For charges separated by a distance L, we have 

S(x) = 2ir5 , (x 2 )i 1 (^+x^r 1 (^-x^ , (4.12) 
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where rj is a step function. When L is much larger than 1/gm the solution of ( [4.1 1| ) can be 
found in the region — L/2 <C X\ <C L/2. In this region ( is approximately independent of X\ 
and satisfies the one- dimensional sine-Gordon equation in x 2 , 

s m 2 sinC = S . (4.13) 



2 



dx^ 

The singular source 5* creates a discontinuity of 27T at x 2 = 0. The solution for x 2 > is 
hence half a sine-Gordon soliton, with £ = 7r at x 2 = 0+ and ( — > as x 2 — > oo. For x 2 < 
it is half the antisoliton with £ = — tt at a; 2 = 0— and ( — > as a; 2 — ► — oo. At distances 
greater than 1/m from the axis, £ vanishes exponentially. Therefore we have ( — > #p 
exponentially. The vortex density m(a;) according to (|4.8| ) vanishes exponentially. Referring 
to ( |4.3| ) we see that the total electric field vanishes exponentially outside the region of width 
1/m, indeed forming a flux tube of thickness 1/m. This is perfectly in accord with our 
expectation for the field profile in a confining theory. 



In ( |4.9|) we see that the energy density of the flux tube is proportional to the energy of 

2 

the sine-Gordon soliton solution with the proportionality coefficient 4-?. The string tension 
therefore is 

2 

a = -^g 2 m . (4.14) 

This is again consistent with the result of the calculation of the spacelike Wilson loop in 
Section 2. We do get a slightly lower value for the string tension (a factor of 2/ir 2 rather 
than 1/4) here than we do from the spatial Wilson loop. This is a consequence of minimizing 
the energy independently in the charged sector. 



C. The massless phase 



The crucial feature of the massless phase is that now D~ l vanishes at zero momentum. 
It is therefore incorrect to expand the classical solution for the dual equation of motion in 
powers of zD. Let us define Xc to be the solution of the classical equation of motion for the 
dual field, 



g 2 D- 1 Xc + zsm( X c + C-dp) = ■ 



(4.15) 



For our purposes it will not be necessary to know the form of \ c explicitly. The essential 
point is that unlike in the massive phase, for the % propagator we must use the complete 
expression 

(xx) connected 



2g 2 D~ L + 2z cos(x c + C - Op) 



(4.16) 



without expanding in z. Calculating the correlation functions of the vorticity via ( 3.15|) , and 
substituting into (|4.4| ), we obtain for the energy 

.2 



AE 



9' 



8vr 2 

r2 



(xc + C)v 2 (x c + C) 



+V 



7T 



2^ 



2a 



Tt D 



-i 



2g 2 D~ L + 2z cos(x c + C - Op) 



-i 



2n-l 



2g'D 



2z 



-l 



(4.17) 

zr 1 . 
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To find the infrared asymptotics we neglect D 1 in the denominators in (|4.17 ) since it 
vanishes at zero momentum. The second term then becomes 

d 2 xV 2 (— • (4.18) 

\cos(xc + C- Op) J 

This is the integral of a total derivative and it vanishes since 1/ cos(x c + C — Qp) does not 
diverge at infinity. 

The profile of ( is therefore determined by minimizing the first term in (|4.17 ). The 
solution is 

Xc + C = 0. (4.19) 

Thus at distances larger than g 2 /z, the induced electric field -^tijdjC, cancels the contribution 
of the vortices ■^-tijdjXc- The total electric field at large distances is therefore just the field of 
the electric dipole. This is again consistent with our expectation for the behavior of electric 
field in a massless phase. 



V. DISCUSSION 

In this paper we studied compact QED3 with mixed action in the gauge invariant Gaus- 
sian variational approach. We find that the theory has two phases. When j3, the coefficient 
of the charge-2 term in the action, is not too negative the theory is massive and confining. 
For f3 large and negative no mass is generated and the theory is not confining. We have also 
calculated variationally the interaction energy of a state with two external charges. 

The authors of [||] studied a minimally modified variational state with no dynamical mass 
within the framework of the pure fundamental theory. As emerges from our analysis the 
vacuum of the pure fundamental theory is massive. This mass has to be introduced explicitly 
in the Gaussian. The nonvanishing string tension is due to this dynamical mass and behaves 
according to the Polyakov formula a ~ g 2 m ~ z 1 ^ 2 where z is the vortex fugacity. 

The massless state is sensible in the theory with mixed action for a > it 2 /A. This 
phase is not confining. The calculation of ||, however, yields a string tension a ~ z which is 
nonvanishing, though parametrically smaller than in the massive phase. We now understand 
that this nonvanishing string tension is an artifact of the minimally modified state. A better 
variational state does not exhibit flux tube formation in the massless phase. 

The main lesson we have learned is that to get a good estimate for the interaction 
energy of external charges we have to introduce a background electric field as an additional 
variational parameter. The minimally modified state of which does not allow for 

this variation is quite poor energetically. It has a wrong profile of the electric field and as 
a consequence overestimates significantly the interaction energy. We expect the same to be 
true in nonabelian gauge theories, so we expect the calculation of the interaction energy of 
external charges there to be more complicated than envisioned in [ |TT| , |T2"f . 

On the positive side we found that the calculation of the spacelike Wilson loop gives a 
good estimate for the interaction potential and the associated string tension even though 
the variational wave function, and indeed the cutoff Hamiltonian theory, is not Lorentz 
invariant. 
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